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Abstract. The MiniBooNE experiment has collected what is currently the world's largest sample of charged current 
single charged pion (CCln^) interactions, roughly 46,000 events. The purity of the CCln^ sample is 87% making this the 
purest event sample observed in the MiniBooNE detector 1 1], The average energy of neutrinos producing CCl?r+ interactions 
in MiniBooNE is about 1 GeV\2], therefore the study of these events can provide insight into both resonant and coherent 
pion production processes. In this talk, we will discuss the long-standing discrepancy in four-momentum transfer observed 
between CClTt^ data and existing predictions. Several attempts to address this problem will be presented. Specifically, the 
Rein-Sehgal j3|,0] model has been extended to include muon mass terms for both resonant |5, 6] and coherent |7] production. 
Using calculations from |8], an updated form for the vector form factor |9] has also been adopted. The results of this 
improved description of CCln^ production will be compared to the high statistics MiniBooNE CCln^ data and several 
existing parametrisations of the axial vector form factor QBIlil] ■ 



INTRODUCTION 

In recent experiments such as K2K lfT2ll , MiniBooNE ifTsIl and SciBooNE ifTill it has been observed that the 
description of the four-momentum transfer in neutrino induced interactions is not sufficient. Currently used Monte 
Carlo generators underestimate the neutrino induced charged current pion production. The predicted shape is different 
from the observed one, especially for the low Q^. There has been an effort in the theoretical community to address 
this discrepancy but mainly focused on data from bubble chamber experiments iHgIITtIi . We will present predictions 
within new models as applied to MiniBooNE Monte Carlo generator 

The K2K collaboration revisited the charged-current coherent pion production, concluding that it is overestimated 
by existing models and they set an upper limit on the coherent contribution to inclusive charged-current interactions 
to 0.6% |il2J . The SciBooNE collaboration [13] set 90% CL upper limits for the same ratio of 0.67% and 1 .36% for 
two samples with mean neutrino energy 1.1 GeV and 2.2 GeV, respectively The MiniBooNE collaboration has 
observed a discrepancy between Monte Carlo prediction and data for the four-momentum transfer in the results for 
CCQE interaction on carbon. In order to fix this problem the MiniBooNE collaboration set an effective value of axial- 
vector mass at M^-^ — 1 .23 ± 0.20Gey, and introduced an effective parameter fc = 1 .019 ± 0.01 1, which modifies the 
Pauli-suppression and reduces the discrepancy between the data and Monte Carlo prediction fis'l. 

The prediction for a neutrino interactions in the MiniBooNE experiment is based on the Monte Carlo generator 
Nuance v3 customized for a carbon nucleus ll25ll . In the Nuance generator the resonance CC\n'^ production is 
describe by the Rein-Sehgal model |3], which is based on the FKR model | T8l[T9l l23l and describes pion production by 
excitation of 18 resonances with a cut at W — 2 GeV in hadronic invariant mass. The resonance production is described 
in terms of helicity amplitudes and the cross section contains interference terms and nonresonant background. The 
vector and axial-vector form factors have a similar form to the ones in case of quasi-elastic scattering, but there 
is an additional factor related to the resonance excitation. Each of the form factors has one free parameter As a 
default for Nuance used in MiniBooNE we use following values: Ma = 1.1 GeV for the axial-vector form factor and 
My = 0.84 GeV for the vector form factor 

The discrepancy in for the CCl7r+ sample is a long standing problem and the main goal of this paper is to 
understand it better. In the first step we will modify the Rein-Sehgal model |3] of resonant pion production to include 
muon mass effects and new form factors. With this new model of pion production we will present predictions for 
CCln^ production in the MiniBooNE detector with various parametrizations of the axial vector form factor For the 
coherent process the effect of mass of final lepton is included by applying the Adler factor lf2l[l . 



In the MiniBooNE detector the struck nucleon is generally below the Cherenkov threshold. That is why the CCln 
sample consists of events produced in following reactions 



V + p - 




+ pH 


- n'^ {resonant production) 


(1) 


V + n - 






- n'^ {resonant production) 


(2) 


v+A - 




+A- 


- 7Z^ {coherent production) 


(3) 



According to the Monte Carlo prediction for the MiniBooNE neutrino beam the biggest contribution comes from 
the excitation of the A(1234) resonance but the noresonant background and coherent production together account for 
about 10% of the sample. In the CCl ;r+ sample 23% of interactions occur on free protons, 46% on bound protons (eq. 
1) and 12% on bound neutrons (eq. 2). Remanning events are coherent production (6%, eq.3) and various background 
interactions (13%). 

In the CCl7r+ reconstruction only one track is reconstructed and assumed to be the jj.^. The muon vertex, track 
angle, and energy, are found with a maximal likelihood fit, with the energy being determined from the total tank 
charge ll22ll . The neutrino energy and four-momentum transfer are reconstructed from the observed muon kinematics, 
treating the interaction as a 2-body collision and assuming that the target nucleon is at rest inside the nucleus and that 
final hadronic state is of the A(1232) resonance. 



FINAL CHARGED LEPTON MASS 



In the Rein-Sehgal models of resonant and coherent pion production via charged current the final lepton is massless. 
Recently, three groups presented extensions of the RS model of resonant production with final charged lepton mass 
included. Kuzmin, Lyubushlin and Naumov |5] (KLN model) showed how one can include lepton mass effects from 
the lepton current. In the Berger and Sehgal model (BRS) |6] in addition to the modification from the KLN model 
the pion-pole contribution is added. We modified our Monte Carlo generator to include the combined KLN and BRS 
models. Predictions for KLN and BRS models are consistent with an independent calculation with the lepton current 
modification and pion-pole contribution by Graczyk and Sobczyk |.23|] . 

For models which introduce a final lepton mass partial cross sections depend on the helicity A 
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KNL showed the form of c^^ ^ for lepton current with the lepton mass. In the limit of massless lepton in the final state, 

(mi 0), the cross section is reduced to the one in FRK model (c^ ' ^ m, cj^ ' ^ v, ' — > 2uv and cj_^' j ^ 0> 61- 

2.18 in ll20ll ) . The helicity amplitudes cr^'^ ^ are calculated as in the FKR model but with three dynamical form factors 
(S, B and C) are modified. 

An axial hadronic current has, in addition to the quark current A^, a pion-pole contribution from the PCAC 
hypothesis, which modifies the axial current. In the case of a massless lepton in the final state the additional term 
vanishes when contracted with lepton current. 

The effect of the lepton mass in the leptonic current is reached by changing only three out of seven dynamical form 
factors (KNL model). Subsequent modification due to the pion-pole term requires alteration of two dynamical form 
factors already modified ^ . 



B^B 
C^C 



^KLN 

(A) 
KLN 
(A) 
KLN 



c 



'BRS 

w 

BRS 



(5) 
(6) 
(7) 



The coherent production in Nuance is described by the Rein-Sehgal model based on by the Adler PCAC theorem 
ll24ll describing neutrino scattering in the forward direction. The extrapolation of the PCAC formula to non-forward 
angles is given by a slowly varying form-factor with w IGeV. 



The important modification of the cross section is due to the final lepton mass. This modification can be found in 
the paper by Adler ll2lll and can be expressed as a simple multiplicative correction factor 

/ 1 0^ \^ I 0^ ■ (0^ -0^ ■ ) 

where 

Qlin = ' eL„ < < 2MEy„a. (9) 

where ymin = m,z/E and ymax = 1 - mi/E. 



FORM FACTORS 

The cross section in the Rein-Sehgal model is expressed in terms of helicity amplitudes but recent analysis on the 
vector and axial forma factors have been presented in terms of form factors in the Rarita-Schwinger model. The 
connection between those two formalism has been provided by Graczyk and Sobczyk |8||], who showed that the vector 
form factor Gy in the RS model can be expressed in term of , C\ and 
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For our analysis we adopted vector form factors obtained by Lalakulich et al. ||^ 

-)2 \ -1 / ra^ -2 



cX 



(-4 



2.13 1 



4M2 



9l 

Ml 



-1.51 1 



4M2 



-1 



Cc = 0.48 1 



4M2 



-1 



Ml 



-2 



0.766M2 



-2 



MW 



(11) 
(12) 

(13) 
(14) 
(15) 



The limit ^ Q gives Gy™(0) — 1.285. This is higher than the value assumed in the RS model. The difference 
comes from the fact that with above shown choice of form factors Cy it is not possible to reproduce the quark model 
prediction of vanishing electric contribution [8]. 

In order to calculate similar relations for the axial vector form factors it was necessary to assume the Adler model 
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Then Ga and Cj are related by | 
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where the coupHng constant C5 (0) is predicted by the PCAC 
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The values of this coupling constant obtained in various quark models are lower than those determined from data. 
The coupling constant evaluation is described in IztIi 

Using the formula ( fTSl l we checked these parametrizations in the Rein-Sehgal model in the Nuance generator. 
In the limit ^ eq. (fTSl l is reduced to the following relation 
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which with W = Ma gives C^(0) = 0.97. 
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FIGURE 1. The axial vector form factor for different parametrizations. . For the form factor in the Rein-Seghal model we 
used the eq. [T8]with W = M^. Left: the parametrizaions which agree better with the ANL data compared to the RS form factor 
Right: form factors which agree better with the BNL data with the RS form factor and dipole C5 with Myj = l.lGeV . Notice that 
for the Lalakulich and Graczyk and Sobczyk vl overlap. 

The Cj (2^) form factor differs in the values of C5 (0), axial mass and the parametrization. Here we present 
predictions using results obtained by three groups which fit the existing data to the axial vector form factor in the 
Rarita-Schwinger formalism. 

Lalakulich et al. : In the results by Lalakulich et al. lIToll the axial form factor has been determined using the model 
of A++(1232) resonance production with contributions only from the amplitude with isospin 3/2. 

In the paper three parametrizations have been presented. For the BNL data the axial form factor from eq.l^Tlhas 
been used (ver 1) with good agreement with the measured distribution for > O.lGeV^. For the ANL data this 
from factor overestimates pion production so two other axial form factors were assumed (eqs.l22landl23Tl and the 
axial mass refitted. 



C^(0) = 1.2, MA^l.OSGeV (21) 
C5(0) = 1.2, MA = l.05GeV (22) 
C^(0) = 1.2, MA = 0.95GeV (23) 



ver.l: C^=C^(0)(^1 
ver.2: Cj^Cj (0)^1 
ver.3: C^=C^(0)(l 




Hernandez et al.: In the model described in ifTlll in addition to the pion production due to A(1232) resonance 
excitation the background terms required by chiral symmetry are included. Assuming the parametrization in 
eq. |24]refitting to the data was necessary since additional background terms changed the normalization. A lower 
value of the axial mass was determined and the correction of the off-diagonal Goldberg-Treiman relation was 
also lower of about 30% (eq.l24li. 

C^=C^(0)(^l+^j (^ + ^) ' <^5 (0)= 0.867, MA=Q.9%5GeV (24) 

Graczyk and Sobczyk: The Rein-Sehgal model has been modified by the authors in order to determine the form 
factor |8] . The fit the dipole form of in the modified Rein-Sehgal model using only the data for scattering 
on proton to avoid the background contribution. They performed two fits for the both sets of data simultaneously 
(ANL and BNL), one for axial mass only (eq. |25] l and one for the axial mass and the Cj (0) coupling constant 
(eq.ESll. 

2 

C^(0) = 1.2, Ma =0.95 GeV (25) 
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C^(0)=0.89, MA = \.\5Gey^ (26) 
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RESULTS 

In figure[T]we present the dependency in the MiniBooNE data for the axial vector form factor Cj parametrizations 
described above and the parametrization from Leitner et al. [28]. For better visualization we divided them into two 
groups, one with better agreement with the ANL data and one with better agreement with the BNL data. Each group is 
compared with the RS form factor obtained from the dipole form using eq. 18, and with dipole with Ma = 1.1 GeV . 
Notice that Cj for the Lalakulich et al. ver. 1 and Graczyk and Sobczyk ver. 1 overlap. The Rein-Sehgal parametrization 
was obtained by applying eq. ( fTsT i to the dipole form of Cj with Ma = 1.1 GeV . 

We have modified the original Rein-Sehgal model in the Nuance generator to include the mass effects and the new 
form factors discussed above (NEW). The resulting predicted distribution are compared to the MiniBooNE data 
in figure |2l We divided prediction into similar groups as in figure [1] In the case of the parametrizations which agee 
with data from ANL we can see good agreement for low < 0.1 GeV^ and underproduction for rest of the region. 
In case of BNL two parametrizations (Lalakulich ver.l and Graczyk and Sobczyk ver. 1) which are almost identical 
but underproduction is not as big as for the ANL results. However, for the dipole form of the with Ma = 1.1 GeV 
agreement for > O.lGeV^ is very good and there is deficit of observed CCl;r+ events for lower Q^. 
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FIGURE 2. MiniBooNE observed CCin^ events as a function of the four-momentum transfer compared with predictions. Top 
figure shows the results with better agreement with the ANL data and bottom figure with BNL data and dipole form of with 
Ma = l.lGeV. For the bottom plot notice that the parametrization by Lalakulich et al. and Graczyk and Sobczyk are virtually the 
same and that the prediction for the dipole form has good agreement with data. 
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